PITHA 98/2 



Magnetic and axial vector form factors as probes of 
orbital angular momentum in the proton 

M. Casu* and L. M. Sehgal* 

00' 

• Institut fur Theoretische Physik (E), RWTH Aachen 

0\ ' 



43 

Oh! 



X 



D-52056 Aachen, Germany 



Abstract 



We have recently examined the static properties of the baryon octet (mag- 
£T) • netic moments and axial vector coupling constants) in a generalized quark 

m ; 

■ model in which the angular momentum of a polarized nucleon is partly spin 

ON 



(S z ) and partly orbital (L z ). The orbital momentum was represented by the 
rotation of a flux-tube connecting the three constituent quarks. The best fit 



is obtained with (S z ) = 0.08 ± 0.15, (L z ) = 0.42 ± 0.14. We now consider 



the consequences of this idea for the g 2 -dependence of the magnetic and ax- 
ial vector form factors. It is found that the isovector magnetic form factor 
Gff vec (q 2 ) differs in shape from the axial form factor Fa{o 2 ) by an amount 
that depends on the spatial distribution of orbital angular momentum. The 
model of a rigidly rotating flux-tube leads to a relation between the magnetic, 
axial vector and matter radii, (r 2 ) mag = / spin (r 2 ) axia i + |/ rb(r 2 )matt, where 
iorb/jspin = \{L x )/Ga, / S pi n + / or b = 1- The shape of F A (q 2 ) is found to be 
close to a dipole with M A = 0.92 ± 0.06 GeV. 
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I. INTRODUCTION 



In a recent paper Q we performed a fit to the magnetic moments of the baryon octet in a 
model in which these quantities are determined partly by the quark spins Am, Ad, As, and 
partly by an orbital angular momentum (L z ), shared between the constituent quarks. The 
model is exemplified by the following ansatz for the proton and neutron magnetic moments: 



fi p = fi u Au + fid Ad + fJ, s As + 
— Hu^d + HdAu + n s As + 



2 1 ' 



2^ + 



(L z ), 

(1) 

(L z ). 



The part containing Am, Ad, As, is the "spin" contribution to the magnetic moments, 
arising from the polarization of quarks and antiquarks in a polarized proton: 

Aq=(q + -q_) + (q + -q_). (2) 

The part proportional to (L z ) is the "orbital" (or convective) contribution, determined by 
the prescription of dividing the orbital angular momentum in proportion to the constituent 
quark masses. The complete set of baryon magnetic moments obtained by this prescription 
is shown in table |I[ These expressions, without the orbital part, were written down in Refs. 

There are two essential elements that go into the above equations for the magnetic 
moments: 

(1) It is assumed that one may use the quark spins Aq in place of the quantities 

8q = (q+ ~ q-) ~ (<?+ - <?-) (3) 

that are appropriate to an expression for the magnetic moment. This approximation 
is justified if antiquarks in a proton carry little polarization. An example of such a 
situation is the chiral quark model ||], in which antiquarks are embedded in a cloud 
of spin-zero mesons. 
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(2) The partition of (L z ) in proportion to the masses of the constituent quarks is based 
on the picture of a baryon as a symmetric three-pronged flux-tube of equal segments 
(figfj), rotating collectively around the spin axis [I]]. The appearence of the same 
magnetons fi u , fid, «s, in the orbital as in the spin part means, in particular, that the 
orbital g-factor has been taken to be g\ = 1. (In a more general description, one could 
interpret (L z ) as (giL z ).) 

II. S z AND L z FROM STATIC PROPERTIES 

The fit to the empirical values of the magnetic moments in |T| was carried out under the 
following constraints: 

(i) The quark magnetic moments were assumed to satisfy fi u = — 2«d, = 0.6/!^. 

(ii) The quark spins Au, Ad and As were constrained to satisfy the measured values of 
the axial vector couplings and a^: 

a (3) = Am - Ad = 1.26 

(4) 

a (8) = Am + Ad - 2 As = 0.58 
These conditions are equivalent to the statement F = 0.46, D = 0.80, in terms of 
which a*- 3 -* = F + D (= Ga, the axial vector coupling constant of neutron decay) and 
a (8) = 3F-D. 

(hi) Each magnetic moment was assigned a theoretical uncertainty of ±0.1ujv (as in Ref. 
H). This ensured that all of the baryons were given essentially the same weight in the 
fit and the x 2 P er degree of freedom was about unity. 

In this manner, the magnetic moments are reduced to functions of three variables, which we 
choose to be fi u , (L z ) and (S z ), the last being defined as 

(S z ) = ^(Au + Ad + As) = ^AJ:. (5) 

The result of the fit is 
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H U = 2.16 ± 0.08, (S z ) = 0.076 ±0.13, (L z ) = 0.42 ± 0.10 (6) 

with x 2 /DOF = 1.1. For the central value of the allowed domain of (S z ) and (L z ) is 
given by the ellipse shown in fig.|2]. Allowing fi u to vary over the interval given in eq.@, we 
obtain the domain shown in fig.|3], from which we infer a final estimate 

(S z ) = 0.08 ± 0.15, (L z ) = 0.42 ± 0.14. (7) 

It is remarkable that the domain of (S z ) and (L z ) determined by the static properties 
of the baryons fulfils rather closely the condition (S z ) + (L z ) = ~. That is, the spin and 
orbital momenta of the quarks and antiquarks saturate the angular momentum of the proton, 
without imposing this as an external requirement. This may be regarded as a posteriori 
justification for the assumption gi = 1. The fact that (S z ) + (L z ) ~ | supports the idea, 
that the spin and orbital angular momentum are linked together by a transition of the form 
q + — > q'_ + M (L = 1), M being a spin-zero meson j|. This in turn provides support to the 
assumption 5q ~ Aq, based on negligible antiquark polarization. 

Also indicated in fig.|3| is the location of two "sign-posts" , that serve as reference points 
in the angular momentum structure: 

(i) NQM: This is the "naive quark model", which describes the nucleon as 3 independent 
quarks in Is orbits, corresponding to (S z ) = |, (L z ) = 0. The SU(6) symmetry of 
the model leads to the prediction Au = |, Ad = — |, As = 0, axial vector couplings 
a^ 3 ' = |, = 1, and the magnetic moment ratio [i v j 'fx n = — |. 

(ii) QPM (As = 0): This is the special case of the quark parton model discussed in 
Ref. 0, in which Au and Ad were allowed to be free, but As was neglected. The 
characteristic prediction of this model is = a^°\ where = Au + Ad + As, 
implying (S z ) = = |(3F — D) = 0.29, the remaining angular momentum being 
attributed to (L z ) = | — (S z ) = 0.21. This version of the QPM leads to the Ellis- Jaffe 
sum rules f| for polarized structure functions: 

/ &{x)dx = \ (F - l -D} , J g?(x)dx = i (F - . (8) 
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In what follows, we consider a test for the presence of orbital angular momentum (L z ) 
and its specific association with the collective rotation of the constituent quarks. 



III. TESTS FOR Lz IN MAGNETIC AND AXIAL VECTOR FORM FACTORS 

We focus on the isovector magnetic moment of the nucleon, obtained by taking the 
difference of \i v and [i n in eq.([l|): 



l^p f^n 

(fi u — fid) 



Ga + \{L z ) 



(9) 



Note that the terms containing As cancel in the difference. We regard this equation as a 
decomposition of the isovector magnetic moment into a part depending on the axial vector 
charge and a part depending on orbital angular momentum. Introducing the abbreviation 

1' _ l^u ~ l^d n f _ f^u ~ 1 ; j v /, n N 

/spin = CjA, Jorb = ^\^z), (,-WJ 

eq.(||) amounts to 

1 = /spin + /orb- (H) 

Returning to the three-parameter fit given by eq.(|6|), we can regard the fitted parameters 
as being (S z ), / spin and / or b (in place of (S z ), fi u and (L z )). For the central value of {S z ), 
the domain of / sp i n and / or b determined by the various magnetic moments is shown in 
fig. |j. The fitted values, taking into account the spread of (S z ), are / sp i n = 0.87 ± 0.03, 
/orb = 0.096 ± 0.03. Considering that these values nearly satisfy the isovector magnetic 
moment relation, / spin + / orb = 1, we use the following approximate values, which fulfil 
eq . (|TTD exactly 

/ spin = 0.90 ± 0.03, / orb = 0.10 ± 0.03. (12) 
Note that the ratio /orb//s P in = Hf^ = | implies (L z ) = 0.42, as given in eq.(0). Eq.(|12|) 
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amounts to the statement, that the isovector magnetic moment is 90% due to quark spin 
polarization and 10% due to quark rotation. 



We now define spatial distributions p maug (r), p ax iai(^) and p rb(r) whose volume integrals 
yield the quantities (p p — p n ), Ga and (L z ) appearing in eq.(^): 



P P - Pn = I d 3 xp mag (r), 

G A = J d 3 X Paxial(r), (13) 

(L z ) = / d 3 xp olh (r). 



The local form of eq.(0) then reads 

Pmag(0 = (flu - fid) 

Introducing, for convenience, " normalized" densities 



PaxialM + ^Porb(^) 



(14) 



Pmag(^) — Pmag (^) /(pp Pn); 

Paxial(r) = p axia \(r)/G A , (15) 

Porb(^) = p oxh {r)/{L z ), 



eq.flTID assumes the form 



Pmag(^) /spin Paxial(^) ~t~ /orb Porb(^)- (16) 



The functions pi(r) all satisfy / p\(r)d 3 x = 1, so that the integrated form of eq.(16) is simply 
the relation (pT|). Fourier transforming eq. ([TB]) , we get a relation between the isovector 
magnetic, axial vector and "orbital" form factors of the nucleon: 

^mag^W 2 ) = /spin -^axial(Q 2 ) + /orb H orh (Q 2 ) (17) 

where 

H i (Q 2 = Q 2 ) = [ pi(r)e i( % (18) 



with ffi(O) = 1. 

The form factor H™™ C (Q 2 ) is an experimentally measured quantity, related to the mag 
netic (Sachs) form factors of the proton and the neutron by 



rrisovcc / r>2\ G P M (Q 2 ) G^^Q 2 ) , 
#mag \Q ) = — — • (19) 

To the extent that G P M (Q 2 ) and G n M (Q 2 ) are both proportional to (1 + Q 2 /0.71 GeV 2 )~ 2 , we 
have 

^g ec (Q 2 ) = -7 M v = 0.84GeV. (20) 
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The (normalized) axial vector form factor is likewise usually parametrized as a dipole 

-^axial(Q 2 ) = — ~2- (21) 

It is clear from eq.(|i~7D, that the difference between H^^ C (Q 2 ) and -ff a xiai(<3 2 ) is a measure 
of the orbital contribution proportional to / or b: in the limit / or b = 0, / spin = 1, these two 
form factors would be identical and we would have Ma = My. 

The orbital form factor H or ^(Q 2 ) is a calculable feature of our model, which ascribes 
the orbital angular momentum to the rigid rotation of a flux-tube. Assuming matter in the 
proton to be distributed as p mat t oc e~ r ^ a , the density of orbital angular momentum p or b is 
proportional to r 2 e~ r l a . The resulting orbital form factor is 

f p iQx r 2 p-r/a J3™ i _ f)2 2 

u in 2 — n 2 \ — ■> — ^ 

HmhiQ ~ Q) - Jr*e-r/a d S x - {1 + Q2 a 2f ^ 

In particular, the rms radius associated with the orbital form factor is 

<r 2 ) orb = ^. (23) 
cr 

This is to be compared with the rms radius of the matter distribution 



(r 2 )matt = ^j, i-e. (r 2 ) orb = ^(r 2 ) mat f (24) 



Eq.(|T7p thus implies a relation between the mean square radii of the various form factors: 



r )mag /spin (?" ) axial + /orb ( r ) 



orb 



2 5 2 

/spin {f )axial + ^/orb )matt- (25) 



To the extent that the matter radius of the proton is assumed to be the same as the magnetic 
radius, we have the prediction 

{t )axial = 7 )mag- (^) 

/spin 

Using the values / or b = 0.10 ±0.03, / spin = 0.90 ±0.03 obtained from the fits to the magnetic 
moments, and the dipole parametrization given in eqs. (|20|) and (pT]) , the above relation yields 

M A = (1.10 ± 0.07) M v = 0.92 ± 0.06 GeV (27) 

in quite reasonable agreement with the value Ma ~ 1.0 GeV deduced from elastic neutrino- 
nucleon scattering 0. It may be remarked here that measurements of elastic pp and pp 
scattering, when interpreted in a geometrical model || tend to give a matter radius slightly 



larger than the charge radius, namely y (r 2 ) matt « 0.89 fm, as compared to y (r 2 ) charge rs 0.84 
fm. If this difference is taken into account, the prediction for Ma obtained from eq.(f26|) 
increases by about one per cent. 

Finally, we can also obtain from eq.([T7|) a more detailed prediction for the shape of 
the axial vector factor H axia [(Q 2 ), in terms of the empirically known magnetic form factor 
#mag eC = (1 + <2 2 /0.71GeV 2 )- 2 and the calculated orbital form factor H orh {Q 2 ) given in 
eq.(p2|). The result is plotted in fig.||, and is close to a dipole with M A ~ 0.92 GeV. 



IV. CONCLUDING REMARKS 

We presented in Ref. [0] a model of the proton as a collectively rotating system of 
quarks, with an orbital angular momentum determined by the baryon magnetic moments 
and the axial vector couplings to be (L z ) = 0.42 ± 0.14. We have now shown that the same 
assumption of a rigidly rotating structure leads to a difference between the normalized axial 
vector and isovector magnetic form factors, which is dependent on the spatial distribution 
of orbital angular momentum. The model of rigid rotation leads to an axial vector form 
factor which is close to a dipole with Ma = 0.92 ±0.06 GeV. Our model of a rotating matter 
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distribution has some similarity to that discussed by Chou and Yang |J, who proposed a 
test for the velocity profile of a polarized proton in hadronic interactions. 
It is of interest to ask how our results for (S z ) and (L z ), namely 



(S z ) = 0.08 ±0.15, (L z ) = 0.42 ±0.14, 



(28) 



compare with those obtained from other considerations. Our fit indicates a dominance of 
orbital over spin angular momentum. This feature is opposite to that in the non-relativistic 
quark model, 



(Sm) = 1> (l„) = o, (NQM) 



(29) 



and closer to the soliton picture of the proton represented by the Skyrme model JT(| 

(S x ) = 0, (L z ) = 1 (Skyrme). 



(30) 



An interesting version of the soliton model,that interpolates between the NQM and Skyrme 



limits, is the chiral quark soliton picture which predicts 

(XQSM). 



F_ _ 5 
D 9 



(31) 



21 + F/D 

In the limit F/D = 5/9, which is the Skyrme model value, one has (S z ) = 0, while in 
the NQM limit Ga = §, F/D = | one has (S z ) = |. For the measured values F = 0.46, 
D = 0.80, this model yields (S z ) = 0.07, which is very close to the estimate in eq. (|28|) . 

Information about (S z ) has also been derived from the analysis of structure functions 
g v { n measured in polarized deep inelastic scattering [p3,|13|] . The integrals of these structure 
functions can be written as 



J gi n (x,Q 2 )dx 



L a (3) + i a (8) 

12 3 



+ 
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a(°)(Q 2 ) 



(32) 



where C± s and Cf are perturbatively calculable coefficients. The singlet axial coupling 
a^(Q 2 ) differs from AS = Am + Ad + As as a consequence of the gluon anomaly. In the 
Adler-Bardeen factorization scheme, a^°\Q 2 ) is related to AE by 



a(°)(Q 2 ) = AS — n^AGfQ 2 ) (33) 

where AG(Q 2 ) is the net polarization of gluons in a polarized nucleon. A determination of 
AS from the measured quantity a^\Q 2 ) is only possible by invoking a model for the polar- 
ized gluon density, and fitting it to the observed Q 2 -dependence of the structure functions. 



The result of one such fit |L4 



is 



(S z ) = -AS = 0.22 ± 0.045 (polarized structure functions) (34) 
2 

Other analyses ( Ui~3| , [J12] , [jn|) obtain values of (5*2) between 0.1 and 0.3. Within errors, 
the result for (S z ) obtained from high energy experiments is compatible with the result ( |28|) 
obtained from a fit to the static properties. 

It remains to be seen whether a specific test of rotational angular momentum (L z ) and 
its radial distribution can be devised. We have argued that the difference in shapes of the 
axial vector and isovector magnetic form factors is a probe of orbital angular momentum. 
A precise determination of Fa{Q 2 ), which does not presume a dipole behaviour from the 
outset, would be of great interest in this respect. 
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Figure Captions 

Fig.l. Flux tube connecting three constituent quarks, rotating collectively around proton 
spin axis. 

Fig. 2. The fitted domain of (S z ) and (L z ) for the central value of \i u . The dotted lines 
represent the baryon magnetic moments of table |I[ 

Fig.3. Allowed domain of (S z ) and (L z ) for the full interval of \x u = 2.16 ± 0.08. 

Fig.4. The fitted domain of / sp i n and / or b for the central value of (S z ). The dotted lines 
represent the baryon magnetic moments of table |I[ 

Fig. 5. Predicted shape of axial vector form factor with / spin = 0.90, / or b = 0.10, compared to 
a dipole with M A = 0.92 GeV. 



Table Caption 

Table I. Parametrization of magnetic moments in the rotating flux-tube model (model (A) of 
Ref. |1J). The fits are based on A = m d /m s = 0.6, fi u = — 2/z rf , fi s = 0.6/i^ (for these 
values, the orbital contribution to the neutral baryons n, 5° and A vanishes). 
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TABLES 



TABLE I. 



Kp) = /^Au + (i d Ad + fi s As + l^ u + \^ d (L z ) 
H{n) = [i u Ad + [i d Au + fi s As + \ix u + \\x d (L z ) 

= HuAu + n d As + n s Ad + j^xl^u + j^X^s (L z ) 
= u u As + u d An + u s Ad + iHx^d + TT5a^ s ^ 
A*(" _ ) = "«As + fl d Ad + fl s Au + 2^AMd+2TAMs (£z) 

At (A ) = i (An + 4Ad + As) + u d ) + § (2Au - Ad + 2As) u s 



+ 
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FIG. 1. 




FIG. 2. 
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FIG. 3. 
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FIG. 4. 
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